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ABSTRACT 



We extend the formalism of the relativistic thermal Sunyaev-Zel'dovich effect to the system 
moving with a velocity /3 = v/c with respect to the cosmic microwave background radiation. In 
the present formalism, the kinematic Sunyaev-Zel'dovich effect for the cluster of galaxies with a 
peculiar velocity (3 is derived in a straightforward manner by the Lorentz boost of the generalized 
Kompanccts equation. We give an analytic expression for the kinematic Sunyaev-Zel'dovich 
effect which is valid up to 0(/3 2 ) with the power series expansion approximation in terms 
of 9 e = ksTe/mc 2 , where T e and m are the electron temperature and the electron mass, 
respectively. It is found that the relativistic corrections to the kinematic Sunyaev-Zel'dovich 
effect are significant. For a typical electron temperature kgT e = lOkcV, one obtains —8.2% 
and +1.3% corrections from the 0((36 e ) and 0{(39l) contributions, respectively. The 0{f3 2 ) 
correction is extremely small, +0.2% for (3 = 1/300 at /csT e = lOkeV. Therefore it can be 
safely neglected. These relativistic corrections are directly reflected on the determination 
of the peculiar velocity (3 of the cluster of galaxies with the observation of the kinematic 
Sunyaev-Zel'dovich effect. 



Subject headings: cosmology: theory — Hubble constant — cosmic microwave background 
radiation — galaxies: clusters: thermal and kinematic Sunyaev-Zel'dovich effects — plasmas: 
Compton scattering 
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1. INTRODUCTION 

Compton scattering of the cosmic microwave background radiation (CMBR) by hot intracluster gas — 
the thermal Sunyaev-Zel'dovich effect — (Zel'dovich & Sunyaev 1969; Sunyaev & Zel'dovich 1972; Sunyaev 
& Zel'dovich 1980a) provides a useful method to measure the Hubble constant H (Gunn 1978; Silk & 
White 1978; Birkinshaw 1979; Cavaliere, Danese, & De Zotti 1979; Birkinshaw, Hughes, & Arnaud 1991; 
Birkinshaw & Hughes 1994; Myers et al. 1995; Herbig et al. 1995; Jones 1995; Markevitch et al. 1996; 
Holzapfel et al. 1997a; Holzapfel et al. 1997b). The Sunyaev-Zel'dovich formula has been derived from 
a kinetic equation for the photon distribution function taking into account the Compton scattering by 
electrons — the Kompaneets equation — (Kompaneets 1957; Weymann 1965). The Kompanects equation 
has been derived with the non-relativistic approximation for the electron. However, the electrons in the 
clusters of galaxies are extremely hot, fceT e = 5 ~ 15kcV (Arnaud et al. 1994; Markevitch et al. 1994; 
Markevitch et al. 1996; Holzapfel et al. 1997a; Holzapfel et al. 1997b). 

Recently attempts have been made to include the relativistic corrections in the Sunyaev-Zel'dovich 
effect (Rephaeli 1995; Rephaeli & Yankovitch 1997). However, it appears that the calculations have not 
been carried out in a manifestly covariant form. For example, equation (4) in Rephaeli (1995), which 
comes from Chandrasekhar (1950), is a non-relativistic formula. Since the extension of the Kompaneets 
equation to the relativistic regime is extremely important in view of many recent measurements of the 
Hubble constant H with the use of the Sunyaev-Zel'dovich effect, we will solve the kinetic equation for the 
photon distribution function in a manifestly covariant form taking into account the Compton scattering by 
electrons. 

Very recently a generalized Kompaneets equation has been derived by several groups: Stebbins (1997); 
Challinor & Lascnby (1998); Itoh, Kohyama & Nozawa (1998). By using the generalized Kompaneets 
equation, analytic expressions for the Sunyaev-Zel'dovich effect have been derived as a power series of 9 e = 
kBT e /mc 2 7 where T e and m are the electron temperature and the electron mass, respectively. It has been 
shown that the results obtained by the power series expansion agree with the previous numerical calculations 
by Rephaeli (1995) and Rephaeli & Yankovitch (1997), thereby proving the validity of their method. It 
is really impressive to see their result agreeing with the result of the manifestly covariant calculation. In 
particular, the convergence of the power series expansion has been carefully studied in the paper by Itoh, 
Kohyama & Nozawa (1998), where the analytic expressions up to 0(6^) have been derived and the results 
have been compared with those of the direct numerical integration of the Boltzmann equation. It has been 
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shown that the power series expansion approximation is sufficiently accurate for the region ksT e < 15kcV 
by taking into account up to O{0\) contributions. 

On the other hand, Sunyaev & Zel'dovich (1980b) have predicted that the CMBR spectrum is also 
distorted if the cluster of galaxies has a peculiar velocity with respect to the CMBR. This effect is called 
the kinematic Sunyaev-Zel'dovich effect. Rephaeli & Lahav (1991) and Hachnclt & Tcgmark (1996) have 
discussed the possibility of the determination of the magnitude of the peculiar velocity by measuring the 
kinematic Sunyaev-Zel'dovich effect. Observation of the kinematic Sunyaev-Zel'dovich effect has been 
done by Holzapfel et al. (1997c) and the limits on the peculiar velocities have been given for the clusters 
Abell 2163 and Abell 1689. At the present stage the error bars are too large to give definite magnitudes 
of the peculiar velocities with the kinematic Sunyaev-Zel'dovich effect. However, more accurate future 
observations are promising in the determination of the peculiar velocities. 

Therefore, relativistic corrections should be carefully taken into account in order to make the 
theoretical expressions for the kinematic Sunyaev-Zel'dovich effect as accurate as that for the thermal 
Sunyaev-Zel'dovich effect. The purpose of the present paper is to calculate the relativistic correction to the 
kinematic Sunyaev-Zel'dovich effect. In the present case, there exist two kinds of the relativistic corrections: 
one is due to the relativisitc electron temperature parameterized by 8 e , and the other is due to the relativistic 
peculiar velocity parameterized by (3 = v jc. In particular we are interested in the interference contributions 
such as O(j30 e ), 0((39 2 ) and 0(f3 2 ). In the present paper we will derive the analytic expressions for the 
thermal and kinematic Sunyaev-Zel'dovich effects from a generalized Kompaneets equation by applying the 
Lorentz boost to the direction of the peculiar velocity of the cluster of galaxies. 

The present paper is organized as follows. In § 2, the general formalism will be presented in deriving 
the thermal and kinematic Sunyaev-Zel'dovich effects on the same footing. With the power series expansion 
approximation, an analytic expression including the relativistic effects will be derived for the thermal and 
kinematic Sunyaev-Zel'dovich effects. Numerical results will be presented in § 3. Finally concluding remarks 
will be given in § 4. 

2. LORENTZ BOOSTED KOMPANEETS EQUATION 

In the present section we will extend the Kompaneets equation to a system (the cluster of galaxies) 
moving with a peculiar velocity with respect to the CMBR. We will formulate the kinetic equation for the 
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photon distribution function using a relativistically covariant formalism (Berestetskii, Lifshitz, & Pitaevskii 
1982; Buchler & Yueh 1976). As a reference system, we choose the system which is fixed to the cosmic 
microwave background radiation (CMBR). The z-axis is fixed to a line connecting the observer and the 
center of mass of the cluster of galaxies (CG). (We assume that the observer is fixed to the CMBR frame.) 
We fix the positive direction of the z-axis as the direction of the propagation of a photon from the cluster to 
the observer. In this reference system, the center of mass of the cluster of galaxies is moving with a peculiar 
velocity (3(= v/c) with respect to the CMBR. For simplicity, we choose the direction of the velocity in the 
x-z plane, i.e. (3 = (p x ,0,/3 z ). 

In the CMBR frame, the time evolution of the photon distribution function n(uo) is written as follows: 



~ = ~ 2 / Wf d3p ' d3k ' w {n{uj)[1 + n ( u 'MW - n ^ + n ^W^ • (2- 1 ) 

2ulu'EE> ' 1 ■ ' 

k = -2(p- k) = -2uE (l - ^-cosaj , (2.4) 

k' = 2(p • k') = 2w'E ( 1 - i^cosa' ) . (2.5) 



In the above W is the transition probability corresponding to the Compton scattering. The four-momenta 
of the initial electron and photon are p = {E,p) and k = (lo, 0, 0, k), respectively. The four- momenta of 
the final electron and photon are p' = (E' ,p!) and k' = (lo' , k'), respectively. The angles a and a' are the 
angles between p and k, and between pand k' , respectively. Throughout this paper, we use the natural unit 
h = c = 1 unit, unless otherwise stated explicitly. Here we note that eq. (2.1) assumes an isotropic photon 
distribution. The kinematic effects lead to anisotropy but to a good approximation eq. (2.1) may be applied 
all the way through the cluster since the levels of anisotropy induced by the kinematic effect are small. 

The electron distribution functions in the initial and final states are Fermi-like in the CG frame. They 
arc related to the electron distribution functions in the CMBR frame as follows (Landau & Lifshitz 1975): 



f(E) = fc(Ec), (2.6) 

f(E') = f c (E' c ), (2.7) 

Ec = 7 (E-p-p) , (2.8) 

E'c = 7 (tf-fi-f) , (2.9) 
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7 = ■ 1 (2.10) 

where the suffix C denotes the CG frame. We caution the reader that the electron distribution function 
is anitotropic in the CMBR frame for (3 ^ 0. By ignoring the degeneracy effects, we have the relativistic 
Maxwellian distribution for electrons with temperature T e as follows: 

Ic{Ec) = e {( E c-m)-(nc-m)}/k B T e + ^ 1 

~ e -{(E c -m)-(nc-m)}/k B T e ^ (2.11) 

where (/zc — m) is the non-relativistic chemical potential of the electron measured in the CG frame. We 
now introduce the quantities 



Ax 



bj' — IV 



k B T e ' 



(2.12) 
(2.13) 



Substituting eqs. (2.6) - (2.13) into eq. (2.1), we obtain 
dn(uj) n f d 3 p 



at 



-2 J AlLd 3 p'd 3 k'Wfc(Ec) [{l + n( W ')}»H 



- { 1 + n(w) } n(u/) e 



i\ Axj(l-0-k') x 7 f3-(k-k') 



(2.14) 



where k and k' are the unit vectors in the directions of k and k', respectively. Eq. (2.14) is our basic 
equation. 

We now expand eq. (2.14) in powers of Ax by assuming Ax <C 1. We obtain the Fokker-Planck 
expansion 



at 



= 2 

+ 2 

+ 2 

+ • • • 

+ 2n 



— h t0 + n(l + n) hs 

— hja + 2(1 + n)— I 2 ,i + n(l + n) ha 
^3 J 3l0 + 3(1 + n)-^ J 3 ,i + 3(1 + n)—h, 2 + n(l + n) J 3 , 3 



where 



(1 + n)Jo + _ Jl + _ j2 + _ j3 + 



J M = ±J^d 3 p'd 3 k'Wfc(E c )(Ax) k e^< k - k 'h e (2.16) 



(2.15) 



Jfc = if / (03^Vd 3 fc'^/c(i?c)(Ax) fc (l - e ^<*- £ '>) . 



(2.17) 
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Analytic integration of eqs. (2.16) and (2.17) can be done with the power series expansion approximation 
of the integrand in terms of the electron momentum p. In Itoh, Kohyama & Nozawa (1998), the systematic 
analysis has been done in order to examine the accuracy of the power series expansion approximation. It 
has been found that the power series expansion approximation is sufficiently accurate for ksTe < 15keV by 
taking into account 0{6^) corrections. 

In addition to 6 e , there is another parameter (3 in eqs. (2.16) and (2.17). For most of the cluster of 
galaxies, (3 -C 1 is realized. For example, (3 « 1/300 for a typical value of the peculiar velocity w=l,000km/s. 
Therefore it should be sufficient to expand eqs. (2.16) and (2.17) in powers of (3 and to retain up to 0(3 2 ) 
contributions. We assume the initial photon distribution of the CMBR to be Planckian with a temperature 

n (X) = ^—, (2.18) 



where 



Substituting the results of eqs. (2.16) and (2.17) into eq. (2.15) and assuming T /T e <C 1, one obtains 
the following expression for the fractional distortion of the photon spectrum: 



An{X) 



+ 



yXe 



x 



,x 



-e e [ Yo + o e Y, + elv 2 + 9% + e 4 e Y 4 ] 



yXe } 



■a 2 



¥£^7 0Pl{0z) [ 1 + OeC! + e 2 e c 2 ' 
— 1 



yXe 



X 



+ V^^(A)[i?o + ^i] 

e A — 1 



(2.20) 



8 = ^ 



cos^ 



Pi(Pz) = B z , 

P*0z) = \ (34 2 - l) , 



(2.21) 
(2.22) 
(2.23) 



where # 7 is the angle between the directions of the peculiar velocity of the cluster {(3) and the initial 
photon momentum (fc) which is chosen as the positive z-direction. The reader should remark that this sign 
convention for the positive ^-direction is opposite to the ordinary one. Thus a cluster moving away from 
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the observer has $ z < 0. However, because of the positive sign in front of the Pi(/3 Z ) term, one obtains 
An(X) < in this case, as one should. The coefficients arc defined as follows: 



Y = -A+X, 

Y 1 = -IQ+^-X-^-X' + ^X* + S 2 



+ 
+ 

+ 
+ 



'_21 7 > 
v 5 5 , 



15 1023 - 



868^ 2 



5 5 



15 2505- 7098~ 2 14253 ~ 3 18594. 



30* 



= 
+ 

+ s : 
s 4 ( 



12059~ 5 128-e 16 - 7 

TuT TT 



140 



7098 



35 

135 



21 " T 105" 

14253- 102267~ 2 156767 1216 - 4 64 - 5 N 

KT + ~ * - -35-* + ^4CT X - — X + Y X \ 

18594 205003 - 1920 ~ 2 1024 -,\ 

h X X 2 H X 3 

^ o«n 7 35 y 



' + 280 
992^ 



a 

1 105 J ' 

30375- 62391-2 614727 -, 124389 - 4 

X X 2 H X A X 4 

128 10 40 10 

-„ 16568 ~„ 7516 ~„ 22 - „ 11 ~„ 



32 ' 128 " 10 " ' 
355703 ^5 16568 - 6 7516 ^ 

80 21 105 7 210" 

- 2 / 62391 614727 - 1368279 - 2 4624139 - 3 157396 - 4 

o -4- X — Jv -4- JC — JC 

' 20 20 80 7 



4U 

7516^.22^8 + ^ 
ins - 210 



+ 



20 

30064 - 5 
— = — X 



2717 - 6 2761 



S 4 - 



7 / ziu / 

124389 6046951 - 248520 - , 

1 x X 2 

160 7 



10 



140 / 

- fi / 70414 465992 - 11792 - , 19778 - ^ 

+ S* 1 X X 2 H X 3 

V 21 105 7 105 , 



481024 -3 15972 ~ 4 
35 7 



682 , 7601 - ^ 

■ A. 



210 



(2.24) 
(2.25) 

(2.26) 



(2.27) 



(2.28) 



C 2 



= 10 



5 



Ix 2 

5 



^S 2 



1117 847^ 2 _ 183^3 + 11^ 
10 10 10 10 

- 2 /847 183 - 121 - 2 \ 11 - 4 



(2.29) 



(2.30) 
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and 



D, = ~l + f Q X, (2.31) 

_ / 19 ~ \ 

£>i = -4+l2X-6X 2 + —X 3 + S 2 (-3 + —X), (2.32) 



y = ot J dlN e , (2.33) 
X = Xcoth [ ^ ] , (2.34) 



sinh I — 



(2.35) 



2 

where N e is the electron number density in the CG frame, the integral in eq. (2.33) is over the photon path 
length in the cluster, and ot is the Thomson cross section. 

Eq. (2.20) is one of our main results. The first line in eq. (2.20) is the expression of the thermal 
Sunyaev-Zel'dovich effect. The relativistic corrections up to 0{9 3 ) was first derived by Challinor & Lasenby 
(1998) and also derived by Itoh, Kohyama & Nozawa (1998) for up to 0(9^) corrections. The second line 
corresponds to the 0{[3 2 ) corrections to the thermal Sunyaev-Zel'dovich effect. The third and fourth lines 
correspond to the kinematic Sunyaev-Zel'dovich effects of the first and second orders in j3, respectively. 
One notes that the first term in the third line of eq. (2.20) corresponds to the lowest-order kinematic 
Sunyaev-Zel'dovich effect (Sunyaev & Zcl'dovich 1980b). The following remarks should be emphasized. 
In the present approach the kinematic Sunyaev-Zel'dovich effects have been derived in a straightforward 
manner from the standard equation for the thermal Sunyaev-Zel'dovich effect by applying the Lorentz boost 
to the initial and final electrons. 

Concerning the photon number conservation of eq. (2.20), the following remarks should be also noted. 
As for the thermal Sunyaev-Zel'dovich effect terms, the photon number conservation is satisfied as discussed 
in Itoh, Kohyama & Nozawa (1998) as well as in Challinor & Lasenby (1998). Namely, Yq, Yi, Yi, Yj 
and Yi terms separately vanish by the integration J dXX 2 . Therefore the first and second lines vanish 
by the integration. On the other hand, the third and fourth lines are propotional to P\{j3 z ) and P2{/3 Z ), 
respectively. Therefore the third and fourth lines vanish by the integration over the solid angle J dfl 7 . 
Therefore the photon number conservation is guaranteed for eq. (2.20). 

Finally we define the distortion of the spectral intensity as follows: 

M = ^^m. (2 .36) 
e x - 1 n (X) 



3. NUMERICAL ANALYSES 



We now study the relativistic contribution for the kinematic Sunyaev-Zel'dovich effect. In Fig. 1 
we have plotted the distortion of the spectral intensity Al/y as a function of X for a parameter set, 
fcsT e =10keV, (3 Z = 1/300. The thermal and kinematic Sunyaev-Zel'dovich effects have been plotted. For 
the kinematic Sunyaev-Zel'dovich effect, the absolute value of the effect is shown. For a receding cluster 
one has AI < 0, whereas for an approaching cluster one has AI > 0. We have shown the decomposition of 
the absolute value of the kinematic Sunyaev-Zel'dovich contribution in Fig. 2. The dashed curve includes 
only the leading order 0(/3) term. The dotdashed curve includes up to the 0(/38 e ) terms. The solid curve 
includes up to O([30l) terms. The dotted curve is the contribution of the 0((3 2 ) terms. It is clear from Fig. 2 
that the relativistic correction is extremely important, which is roughly 8% of the leading order contribution 
of the kinematic Sunyaev-Zel'dovich effect. It is also clear from Fig. 2 that the 0((3 2 ) correction is very 
small and it can be safely neglected. 



3.1. RAYLEIGH JEANS REGION 



In the Rayleigh- Jeans limit where X — > 0, eq. (2.20) is further simplified: 



An{X) 

~MX) 



-2yO e 
-2yP 2 
+ ypPi(P z ) 



6 e 



123 
"40" 



1989 3 
280 e 



14403 f 
640 



10 



1 



2 
5 

15 5 



13 



(3.1) 



For this expression also, one should note the sign convention for the positive z-axis stated below eq. (2.23). 
For the purpose of illustration we have calculated the kinematic Sunyaev-Zel'dovich effect for a parameter 
set ksTe = lOkeV, (3 Z — 1/300. In this region, the leading order kinematic Sunyaev-Zel'dovich effect is 
8.8% of the thermal Sunyaev-Zel'dovich effect. The relativistic corrections of the O((30 e ) and 0(/3 2 ) terms 
arc 0.07% and 0.02% of the thermal Sunyaev-Zel'dovich effect, respectively. Therefore the relativistic 
corrections of O(/30 e ) and 0(/3 2 ) are safely neglected. 
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3.2. CROSSOVER FREQUENCY REGION 



It is well known that the kinematic Sunyaev-Zel'dovich effect becomes extremely important in the 
crossover frequency region, where the thermal Sunyaev-Zel'dovich effect vanishes. Therefore the accurate 
position of the crossover frequency X is extremely important. In Itoh, Kohyama & Nozawa (1998), the 
following fitting formula for X has been obtained. 



The errors of this fitting function are less than 1 x 10 3 for < k B T e < 50keV. 

Substituting eq. (3.2) into eq. (2.36), we have calculated the kinematic Sunyaev-Zel'dovich effect at 
X = Xq as a function of ksT e . The result has been plotted in Fig. 3. The dashed curve is the leading 
order 0((3) contribution. The dotdashed curve includes up to the 0{[36 e ) contribution. The solid curve 
includes up to the 0(@6 2 ) contribution. The 0{f3 2 ) contribution is very small and is again safely neglected. 
It is clear that the relativistic corrections become significant as the electron temperature increases. At 
a typical electron temperature fcgT e = lOkeV, one obtains the —8.2% and +1.3% corrections from the 
0{[39 e ) and 0(f36 2 ) contributions, respectively. The 0((3 2 ) correction is extremely small. It is +0.2% 
for (5 = 1/300 at fcsT e = lOkeV. Therefore one can safely neglect the 0{[3 2 ) corrections. Although the 
thermal Sunyaev-Zel'dovich effect is exactly zero at X = Xq, the slope is very large as shown in Fig. 1. 
At ksTf, — lOkeV, for example, the thermal Sunyaev-Zel'dovich effect becomes ±17% of the kinematic 
Sunyaev-Zel'dovich effect if one shifts X by ±1%. Therefore the precise determination of the position of 
Xo is extremely crucial in isolating the kinematic Sunyaev-Zel'dovich effect. 

It is well known that the kinematic Sunyaev-Zel'dovich effect is useful in determining the peculiar 
velocity; see for example, Sunyaev & Zel'dovich (1980b) and Rephaeli & Lahav (1991). The relativistic 
corrections are directly reflected on the determination of the peculiar velocity of the cluster of galaxies with 
the observation of the kinematic Sunyaev-Zel'dovich effect. The crossover frequency is shifted toward higher 
(lower) frequency region by the kinematic Sunyaev-Zel'dovich contribution when the cluster is moving 
outward (inward) with respect to the observer. The shift is given by AA = X — X ,p, where X is given 
by eq. (3.2) and X ,p is the new crossover frequency in the presence of the kinematic Sunyaev-Zel'dovich 
effect. We define a function h as follows. 



In Fig. 4 we have plotted ft. as a function of fcsT e for the typical peculiar velocities v = 500km/s, l,000km/s 




(3.2) 




(3.3) 
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and l,500km/s, which correspond to (3 = 0.5/300, 1/300 and 1.5/300, respectively. These curves show an 
excellent scaling for 5keV < fcsT e < 15keV. The function h can be fitted as follows. 

h = — ^ + - f 3 + a 3 + a 4 9 e + a 5 8 2 e , (3.4) 

Oe "e.mm (U e — U e ^ m i n ) 

where 8 e ,min is the minimum value of 9 e which has a real solution for a given (3. We fix 6 e ^ min — 1.654x 10~ 3 . 
The other constants are di = 3.857xl0~ 3 , a 2 = -4.631 x 10~ 6 , a 3 = 1.370xl0~ 2 , a 4 = 1.014xl0~ 2 , 
and as = l.OOOxlO -2 . The errors of this fitting function are less than 4% for 5 < ksT e < 15keV and 
0.5/300 < (3 < 1.5/300. 



4. CONCLUDING REMARKS 

By applying the Lorentz boost to the standard formalism of the extended Kompaneets equation, we 
have derived a formalism of the kinematic Sunyaev-Zcl'dovich effect for the cluster of galaxies with a 
peculiar velocity (3. With the power series expansion approximation in terms of the electron temperature 
6 e and the peculiar velocity (3, we have derived an analytic expression for the kinematic Sunyaev-Zel'dovich 
effect which includes the relativistic corrections of up to O((30 2 ) and 0(/3 2 6 e ). It has been found that 
the relativistic correction is significant. At a typical temperature ksT e = lOkcV, one obtains the —8.2% 
and +1.3% corrections from the O((30 e ) and 0((39 2 ) contributions, respectively. The 0([3 2 ) correction is 
small. It is +0.2% for (3 = 1/300 at ksTe = lOkcV. These relativistic corrections are directly reflected on 
the determination of the peculiar velocity of the cluster of galaxies with the observation of the kinematic 
Sunyaev-Zel'dovich effect. Finally we have fitted the shift of the crossover frequency AA+,. The shift has an 
excellent scaling for 5 < k B T e < 15keV and 0.5/300 < (3 < 1.5/300. The errors of the fitting function are 
less than 4%. 

We received a preprint "Cosmic Microwave Background Radiation in the Direction of a Moving Cluster 
of Galaxies with Hot Gas: Relativistic Corrections" from S. Y. Sazonov and R. A. Sunyaev after the 
submittal of our original manuscript to The Astrophysical Journal. Most of our original results agreed with 
theirs (their equation (12)) except for the term proportional to f3 2 (the second term on the right-hand-side 
of equation (2.20) in the present paper). Although this term makes a negligible contribution in the real 
situation of the clusters of galaxies, we endeavored to find out the cause for the discrepancy. Later, this 
discrepancy was pointed out by our referee A. Challinor who carried out an independent calculation. 
Finally, we have found out an error in our original calculation and have corrected it. Now our result agrees 
with that of Sazonov and Sunyaev (although they stop at a lower order of expansion) as well as that of 



-13- 



Challinor. We sincerely wish to thank S. Y. Sazonov, R. A. Sunyaev, and A. Challinor for motivating us to 
obtain the correct term proportional to 1 . We also wish to thank R. A. Sunyaev for his very informative 
communication to us. Last but not least, we wish to thank A. Challinor for his invaluable work as our 
referee, examining our mathematical results thoroughly and giving us many useful suggestions. 
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Figure Captions 

• Fig.l. Spectral intensity distortion Al/y as a function of X for fcsT e = lOkeV, (3 = /3 2 =l/300. 
The thermal and kinematic Sunyaev-Zeldovich effects have been plotted. For the kinematic 
Sunyaev-Zel'dovich effect, the absolute value of the effect is shown. 

• Fig. 2. The decomposition of the absolute value of the kinematic Sunyaev-Zel'dovich contribution for 
k B T e = lOkeV, (3 = /3 z =l/300. The dashed curve is the leading order 0((3) term. The dotdashed 
curve is the 0{(3) + 0{f39 e ) terms. The solid curve is the 0(/3) + O{j30 e ) + 0((36l) terms. The dotted 
curve is the contribution of the 0{f3 2 ) terms, where a factor 10 was multiplied in order to make the 
curve visible in this figure. 

• Fig. 3. The absolute value of the kinematic Sunyaev-Zel'dovich effect at X = X as a function of 
ksTe for (3 Z — 1/300. The dashed curve is the leading order O(fi) term. The dotdashed curve is the 
0(0) + 0{/36 e ) terms. The solid curve is the 0((3) + 0{[36 e ) + 0{f36 2 e ) terms. 

• Fig. 4. The shift h of the crossover frequency as a function of fcgT e . The dashed, dotdashed and solid 
curves correspond to v = 500km/s, l,000km/s and l,500km/s, respectively. 



